arXiv:hep-ph/0211085v4 11 Oct 2005 


FERMILAB-PUB-02/274-T 
ANL-HEP-PR-02-094 


Exclusive Double-Charmonium Production 
from e+e“ Annihilation into a Virtual Photon 

Eric Braaten 

Physics Department, Ohio State University, Columbus, Ohio 43210 and 

Fermi National Accelerator Laboratory, 

P. 0. Box 500, Batavia, Illinois 60510 

Jungil Lee 

Department of Physics, Korea University, Seoul 136-701, Korea 
(Dated: February 1, 2008) 

Abstract 

We calculate the exclusive cross sections for e^e“ annihilation into two charmonium states 
through a virtual photon. Purely electromagnetic contributions are surprisingly large, changing 
the cross sections by as much as 21%. The predicted cross section for J/fj + r?c is about an order 
of magnitude smaller than a recent measurement by the BELLE Collaboration, although part of 
the discrepancy can be attributed to large relativistic corrections. The cross sections for S'-wave + 
P-wave, P-wave + P-wave, and S'-wave + D-wave charmonium states are also calculated. R may 
be possible to discover the D-wave state r/c 2 (lP) at the B factories through the mode J/'f ric 2 , 
whose cross section is predicted to be about a factor of 10 smaller than J/fj P tie- 
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I. INTRODUCTION 


NRQCD factorization is a systematic framework for calculating the inclusive cross sections 
for producing heavy quarkonium [1]. The cross section for a charmonium state is expressed 
as the sum over cc channels of products of perturbative cc cross section and nonperturbative 
NRQCD matrix elements. The relative importance of the various terms in the factorization 
formula is determined by the order in of the cc cross section, kinematic factors in the 
cc cross section, and the scaling of the NRQCD matrix element with the relative velocity 
V of the charm quark. Among those NRQCD matrix elements that scale with the minimal 
power of V is the one associated with the color-singlet cc channel whose angular momentum 
quantum numbers match those of the charmonium state H. The old color-singlet model for 
quarkonium production [ 2 ] can be dehned by keeping this channel only. 

If the charmonium is the only hadron in the initial or hnal state, the color-singlet model 
should be accurate up to corrections that are higher order in v. The simplest examples of 
such processes are electromagnetic annihilation decays, such as J/fj —>• and rjc 77 , 

and exclusive electromagnetic production processes, such as 77 — rjc- Another process 
for which the color-singlet model should be accurate is e’'“e“ annihilation into exactly two 
charmonia. There are no hadrons in the initial state, and the absence of additional hadrons 
in the hnal state can be guaranteed experimentally by the monoenergetic nature of a 2 -body 
hnal state. For many charmonia H, the NRQCD matrix element can be determined from 
the electromagnetic annihilation decay rate of either H or of another state related to H 
by spin symmetry. Cross sections for double-charmonium can therefore be predicted up to 
corrections suppressed by powers of without any unknown phenomenological factors. 

One problem with annihilation into exclusive double charmonium is that the cross 
sections are very small at energies large enough to have conhdence in the predictions of 
perturbative QCD. A naive estimate of the cross section for J/-\- rjc in units of the cross 
section for is 


R[J/^ + 7c] ~ ■ (1) 

The 2 powers of a* are the fewest required to produce a cc-f cc hnal state. There is a factor 
of {rUcvY associated with the wavefunction at the origin for each charmonium. These factors 
in the numerator are compensated by factors of the beam energy A'beam in the denominator 
to get a dimensionless ratio. As an example, consider e^e“ annihilation with center-of-mass 
energy 2i?beam = 10.6 GeV. If we set 0.3, as ~ 0.2, and rUc ~ 1.4 GeV, we get 

the naive estimate R[J/ip + 7 c] ~ 4 x 10“^. This should be compared to the total ratio 
R[hadrons] ~ 3.6 for all hadronic hnal states [3]. The decay of the J/ijj into the easily 
detectable e+e“ or modes suppresses the observable cross section by another order of 

magnitude. 

Fortunately, the era of high-luminosity B factories has made the measurement of such 
small cross sections feasible. The BABAR and BELLE detectors have each collected more 
than 10^ continuum e’''e“ annihilation events and more than 10® events on the T(45'), 75% 
of which are continuum e+e“ annihilation events. The BELLE Gollaboration has recently 
measured the cross section for J/ip -\- rjc [4]. They also saw evidence for J/ip Xco 

and J/Ip + riciplS) events. 

In this paper, we calculate the cross sections for exclusive double-charmonium production 
via annihilation into a virtual photon. This process produces only charmonium states 
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with opposite charge conjugation. The cross sections for charmonium states with the same 
charge conjugation, which proceed through e’''e“ annihilation into two virtual photons, will 
be presented in subsequent papers [5]. We carry out the calculations in the color-singlet 
model including not only the diagrams of order but also the purely electromagnetic 
diagrams of order which are surprisingly large. Our result for the cross section for 
J/ip + rjc is about an order of magnitude smaller than the recent measurement by the 
BELLE Collaboration, although part of the discrepancy can be attributed to large relativistic 
corrections. The cross sections for 5'-wave -|- P-wave, P-wave -|- P-wave, and S'-wave -|- D- 
wave charmonium states are also calculated. The cross section for J/ip + r]c 2 {^D) is predicted 
to be about a factor of 10 smaller than for J/iIj + rjc, which may be large enough for the 
P-wave state ? 7 c 2 (lP) to be discovered at the B factories. 


II. COLOR-SINGLET MODEL CALCULATIONS 


In this section, we use the color-singlet model to calculate the cross sections for 
annihilation through a virtual photon into a double-charmonium hnal state H 1 + H 2 . Charge 
conjugation symmetry requires one of the charmonia to be a C* = — state and the other to 
be a C* = -t- state. The C = — states with narrow widths are the = 1 states J/ijj and 
'ip{2S), the I"*"" state he, and the yet-to-be-discovered 2 state ■^2(1-0). The C = + states 
with narrow widths are the states rjc and rjei^S), the states Xcj(l-P), J = 0,1,2, 
and the yet-to-be-discovered 2“+ state ? 7 c 2 (lP). We express our results in terms of the ratio 
R[Hi + P 2 ] dehned by 


R[Hi + P2] 


(T[e+e ^ Hi + H2 
(T[e+e“ — 


( 2 ) 


In the text, we give only the results for R summed over helicity states. In the Appendix, 
we give also the angular distribution dR/d cos 6 for each of the helicity states of Hi and 
H 2 . These results may facilitate the use of partial wave analysis to resolve the experimental 
double-charmonium signal into contributions from the various charmonium states. 


A. Asymptotic behavior 


When the beam energy Pbeam is much larger than the charm quark mass rric, the 

relative sizes of the various double-charmonium cross sections are governed largely by the 
number of kinematic suppression factors r^, where the variable r is dehned by 


r 


2 


4m^ 

/ 7'2 

■^beam 


( 3 ) 


If we set nic = 1.4 GeV and Pbeam = 5.3 GeV, the value of this small parameter is = 0.28. 

The asymptotic behavior of the ratio R[Hi + H 2 ] as r —0 can be determined from the 
helicity selection rules for exclusive processes in perturbative QCD [6, 7]. For each of the 
cc pairs in the hnal state, there is a suppression factor of due to the large momentum 
transfer required for the c and c to emerge with small relative momentum. Thus, at any 
order in a*, the ratio R[Hi + H 2 ] must decrease at least as fast as as r —0. However 
it may decrease more rapidly depending on the helicity states of the two hadrons. There 
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is of course a constraint on the possible helicities from angular momentum conservation: 
|Ai — A2I = 0 or 1 . The asymptotic behavior of the ratio R[Hi{\i) +-^^2(^2)] depends on the 
helicities Ai and A2. The helicity selection rnles imply that the slowest asymptotic decrease 
R ^ can occnr only if the snm of the helicities of the hadrons is conserved. Since there are 
no hadrons in the initial state, hadron helicity conservation reqnires Ai + A2 = 0 . The only 
helicity state that satisfies both this constraint and the constraint of angular momentum 
conservation is (Ai,A2) = ( 0 , 0 ). For every nnit of helicity by which this rnle is violated, 
there is a further suppression factor of r^. The resulting estimate for the ratio R at leading 
order in ag is 


Rqcb[Hi{Xi) + H2{X2)] ~ (4) 

The factor of for a charmoninm state with orbital angular momentum L comes from the 
NRQCD factors. At leading order of ag, there may of course be further suppression factors 
of that arise from the simple structnre of the leading-order diagrams for e’''e“ —> cci -f cci 
in Fig. 1 , bnt these snppression factors are nnlikely to persist to higher orders in ag. 

The QED diagrams for > cci(^S'i) +cci in Fig. 2 give contribntions to R[J/'ip + H2\ 

that scale in a different way with r. As r — 0 , the contribntion to the cross section from 
these diagrams factors into the cross section for 'j + H2 and the fragmentation fnnction for 
7 — >• J/'ip. This fragmentation process prodnces J/ip in a Xj/^ = ±1 helicity state. The 
hard-scattering part of the process prodnces only one cc pair with small relative momentnm, 
so there is one fewer factor of relative to Eq. ( 4 ). The cross section for 7 -|- FTi is still 
snbject to the helicity selection rules of perturbative QCD, so the pure QED contribution 
to the ratio R has the behavior 

RQED[W(±l)+hf2(A2)] ~ (5) 

There may also be interference terms between the QCD and QED contribntions whose 
scaling behavior is intermediate between Eqs. ( 4 ) and ( 5 ). 


B. Projections onto charmoninm states 

The 4 QCD diagrams for the color-singlet process 7* —cci -|- cci are shown in Fig. 1 . 
We take the npper cc pair in Fig. 1 to form a C = — charmoninm Hi with momentnm Pi 
and the lower cc pair to form & C = + charmoninm H2 with momentum P2. There are also 
QED diagrams for 7* —cci -|- cci that can be obtained from the QCD diagrams in Fig. 1 
by replacing the virtnal glnons by virtnal photons, bnt they are snppressed by a factor of 
a/ag. However if one of the charmonia is a 1 state like a J/t/’, there are the additional 
QED diagrams in Fig. 2 . Althongh they are also snppressed by a factor of a/ag, they are 
enhanced by a kinematical factor of 1/r^ and therefore can be more important than one 
might expect. 

To calculate the matrix element for e+e“ — Hi{Pi) + H2{P2), we start from the matrix 
element for e+e“ —> c{pi)c{pi) + c{p2)c{p2) with the charm qnarks and antiqnarks on their 
mass shells: pj = pj = Tn^. For each of the cc pairs, we express the momenta in the form 

p = \P + q, 
p = \P -q, 
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(6a) 

(6b) 







where P is the total momentum of the pair and g is a relative momentum that satishes 
g ■ P = 0. If the cc pair is in a spin-singlet color-singlet state, the matrix product of the 
Dirac and color spinors for the c and c can be expressed as [9] 

v{p)u{p) = ^ - mP) 75 (f + 2E) {i) + 1 

A^/2E{E + mc) \VNc 

where = P^/4 = — g^, Nc = 3, and the last factor involves the unit color matrix 1 . 

If the cc pair is in a spin-triplet color-singlet state, 75 in Eq. (7) is replaced by where es 
is a spin polarization vector satisfying ■ ej = — 1 and P ■ 65 = 0. 



5 



In the spin-singlet case, the expansion of the matrix element in powers of q has the form 

M [cc{S = 0)] = ^ + + C^rq^q^ + .... (8) 

The matrix elements at leading order in v for the spin-singlet charmonium states rjc, hc{lP), 
and ? 7 c 2 ( 1 -D) can be read off from this expansion; 

(Oi), ^ 


M[ric\ = 

M[K{\)] = 
M[ric2{\)] = 


'rjc 


2Ncmc 

{Oi)k 

2Nrml 


A, 


1/2 




(Oi) 


ric2 


2Nrml 


1/2 


a 


^ J'Cr/i J'TIV _J_ JT^ _ ^JCTTJfliy 




(9a) 

(9b) 

(A), (9c) 


where e'^(A) in Eq. (9b) is a spin-1 polarization vector and e^y{\) in Eq. (9c) is a spin-2 
polarization tensor. The tensor in Eq. (9c) is 

pfi pu 

-r + ( 10 ) 


I 


Ami ’ 


where P is the expansion of p -|- p to leading order in g, so it now satisfies P^ = Am^. The 
NRQCD matrix elements (Oi),,^ and {Oi)h^ in Eqs. (9a) and (9b) are the vacuum-saturated 
analogs of the NRQCD matrix elements (Oi(^S'o))r,c and (Oi(^Pi))/i^ for annihilation decays 
dehned in Ref. [1]. The NRQCD matrix element (Oi),j ^2 in Eq. (9c) is the vacuum-saturated 
analog of an NRQCD matrix element (Oi(^D 2 ))jjc 2 which in the notation of Ref. [1] is dehned 
by 

(OiQD2 ))^.. = (hc2|^^(-f)'B 0 D 0 D ^'Vlhc2). (11) 

A projection onto the D-wave state that is equivalent to Eq. (9c) but expressed in terms of 
wavefunctions at the origin is given in Ref. [10]. 

In the spin-triplet case, the expansion of the matrix element in powers of q has the form 

M[cc{S = 1)] = (Alp + + Cp^rq^q^ + ...) e§. (12) 

The matrix elements at leading order in v for the spin-triplet charmonium states J/i/’, 
Xcj(lP), t/’i(l-D); and i/’ 2 (l-D) can be read off from this expansion: 

(Oi)j/i/.^ 


MWiX)] = 

M[xco] = 

-^[Xci(A)] = 
M[Xc2{X)] = 
-^[^i(^)] = 


2Ncmc 

(Qi)xco 

2Ncml 

(Qi)xcb 

2Ncml 

(Qi)xc 2 

2Ncml 

2Ncml 

{Qi)^2 

2NpiP 


^pC^(A), 


1/2 


B _ 

Vs ’ 


1/2 


B, 


per 


2mcV2 




1/2 




pa 


T 1 

- 1 - _ Jpf^JPi' 

2 3 




(A), 


1/2 


c 


par 


— (P'^P^ -1- jpTj^p'j _ }ljppj 
2 3 


1/2 


C 


par 


2mcVQ 




(13a) 

(13b) 

(13c) 

(13d) 

(13e) 

(13f) 
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The NRQCD matrix elements (Oi) and (Oi)xcj I'he vacuum-saturated analogs of the 
NRQCD matrix elements (Oi(^S'i))j/^ and (Oi(^-Pj))xcj annihilation decays defined in 
Ref. [1]. The NRQCD matrix elements (Oi)^j and (Oi )^2 are the vacuum-saturated analog 
of the NRQCD matrix elements (Oi(^Di))^^ and (Oi(^D 2))^2 which are dehned by 

(OiQD2))^2 = (14b) 


C. S-wave -|- S-wave 


The matrix element for e~^{k 2 )e {ki) J/ijj + r]c can be written as 

2 

M = -^^v{k2)'y'"u{ki) {J/^jJ + r]c\J^,{0)\^), (15) 

where Cc = -|-| is the electric charge of the charm quark and s = is the square of the 

center-of-mass energy. Upon simplifying the vacuum-to-J/'^ -|-?7c matrix element, it reduces 
to the general form required by Lorentz covariance: 

(J/^(Pi,e)+r/e(P2)|JM(O)|0) = tAe^,xaPiP2e% (16) 


where the coefficient A is 

1287rQ;o 


A = 




1/2 /iVf - 1 e^a , 1 e^a 


2m 


+ 


NrU. a 


After squaring the amplitude and integrating over phase space, we obtain 
for the ratio R defined in Eq. (2): 

R[j/^ + x\r‘^ - Yfr\l - (Oi) ^ 

9 mk 


(17) 

our hnal result 

(18) 


where the coefficients X and Y are 

X 

U 



(19a) 

(19b) 


If we set as = 0.21, their numerical values are X = 0.450 and Y = —0.0344. Note that the 
ratio (18) depends on the charm quark mass rric explicitly and also through the variable r 
dehned in Eq. (3). The a'^a'j. term in the cross section for e+e“ — J/'ijj + rjc was calculated 
previously by Brodsky and Ji [8]. They presented their result in the form of a graph of R 
versus 1/r^, but they did not give an analytic expression for the cross section. 

The only helicity states that contribute to Eq. (18) at this order in are (Ai, A 2 ) = 
(±1,0), which violate hadron helicity conservation by 1 unit. The QCD contribution to R 
scales like in accord with Eq. (4). There is no reason to expect the amplitude for 


7 



the hadron-helicity-conserving state (0, 0) to vanish at next-to-leading order in a.^, so the 
asymptotic behavior of the QCD contribution as r —0 is probably R ~ 

The pure QED contribution from the diagrams in Fig. 2 scale like in accord with 

Eq. (5). The interference term is suppressed only by a/{asr‘^)-, so the QED effects are larger 
than one might expect. If we set ^/s = 10.6 GeV and rric = 1.4 GeV, the electromagnetic 
correction increases the cross section by 29%. 

The factor of (1 — in Eq. (18) is the nonrelativistic limit of (TciviZ-Ebeam)^, where 

PcM is the momentum of either charmonium in the center-of-momentum frame. It can be 
expressed as 

PcM ^ 

■^beam "S 

where \{x, y, z) = -[- z'^ — 2{xy + yz + zx). In (18), one factor of PcM/-Ebeam comes 

from the phase space for Hi + H 2 -, while the other two come from the square of the amplitude 
(16). 


D. S-wave -|- P-wave 


The ratios P for J/'ip + Xcj are 



TT^CK^ 

RlJ/i- + Xcj] = 

X^Fj{r, Yyil (Oi)J rWxcj ^ 

mp. 

(21) 

where the functions Fj{r,Y) are 




Fq\ 

XY) 

= 2[4V-6(V + 3)r 

2 + 7r^; 

P + r2[4 + 2(V + 5)r2-3r^]^ 

(22a) 

Pil 

X Y) 

= 3[8Y - 2Yr‘^ + 

]' + 3% 

■[AY + 2iY + 2)r‘^ - 3ry 




+ 3r^[2(3V + 4) - 


(22b) 

P 2 I 

X Y) 

= [8Y-6{Y + 2y 

-|- llr^ 

f + 2ry + 2{Y - - 3ry 




+ 3r2[4V - 2(V + 2y + Sr^ + 3r^[2(V + 2) - br-y + 6r^[2V - 

r2]2. (22c) 


These expressions have been expressed as sums of squares of terms that correspond to the 
helicity amplitudes. For Xco and Xc 2 , there are QGD contributions from all the helicity states 
(Ai, A 2 ) compatible with angular momentum conservation, so the leading contribution scales 
like For Xci, there is no contribution from the hadron-helicity-conserving state (0, 0), 

so the QGD contributions are suppressed by relative to those for Xco and Xc 2 - The pure 
QED contributions to Xcj scale like a^v^r'^ for all J in accordance with Eq. (5). The QED 
contribution is suppressed by a factor of for J = 0 and 2 but only by /{alr^) 

for J = 1. If we set ^/s = 10.6 GeV and me = 1.4 GeV, the QED corrections change the 
cross sections by -|-5.0%, —5.5%, and -|-11% for J = 0, 1, and 2, respectively. 

The ratio R for rje + he is 


Rh + h.] = . 


(23) 


where the function H{r) is 


H{r) = 2r^(r^ — 2)^ -|- (3r^ — 6r^ -|- 4)^. 
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( 24 ) 



The dependence on a appears only in the overall factor because the QED contribution 
conies from diagrams with the same topology as the QCD diagrams in Fig. 1 of Ref. [22]. The 
QED contribution from the photon-fragmentation diagrams in Fig. 2 of Ref. [22] vanishes 
because the -|- parity of he does not allow the direct coupling to a single virtual photon. The 
pure QCD contribution to the result in Eq. (23) is obtained by substituting X = 4/9. This 
result was first calculated correctly in Ref. [21]. 


E. P-wave -|- P-wave 


The ratios R for he + Xcj are 


R[hc + Xcj] — 
where the functions Gj{r) are 


108 

X‘^Gj{r)r\l ^ 2 ) 3/2 

m]X 

(25) 

Goir 

) = 2r^{6-r^f, 

(26a) 

G'i(r 

) = 24±78r2(2-rQ^ 

(26b) 

G' 2 (r 

) = 3r2(4-5r2)^±7r®. 

(26c) 


At this order in there is no contribution to the cross sections for Xco and Xc 2 from the 
helicity state (0, 0), so the ratios R[he + Xcj] for Xco and Xc 2 are suppressed relative to that for 
Xci by a factor of r^. The QED contribution increases the cross section by 2a/(3as) ~ 2%. 


F. S'-wave -|- P-wave 

The ratio R for + ? 7 c 2 is 

am + 1.2} = ^ X^(r - 2 rrr\l - , ( 27 ) 

At this order in Og and a, the only helicity states that contribute are (±1,0). Thus the 
QCD contribution to the ratio R scales like in accord with Eq. (4), while the pure 

QED contribution scales like in accord with Eq. (5). If we set y/s = 10.6 GeV and 

nie = 1.4 GeV, the QED correction increases the cross section by about 15%. 

The ratio R for ± r]e is 

+ xd = ^ X'^(26Y - 21r^ + 10r‘) V^(l - (jg) 

At this order in ag and a, the only helicity states that contribute are (±1, 0). Thus the QCD 
contribution to the ratio R scales like in accord with Eq. (4), while the pure QED 

contribution scales like a‘^v^^r‘^ in accord with Eq. (5). The QED contribution increases the 
cross section by 41%. 

The ratio R for 1^2 ± Vc is 

Rli’2 + »)J = ^ (OdfejOi)’.- , (29) 

54 rrie 

There is a contribution from the helicity state (0, 0) that satishes hadron helicity conserva¬ 
tion, so the ratio R scales like in accord with Eq. (4). The QED contribution increases 
the cross section by 2a/{3ag) ~ 2%. 
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TABLE I: The NRQCD matrix elements {Oi)h for the charmonium states H in units of (GeV)^^'*'^ 
where L = 0,1,2 for S-wave, P-wave, and D-wave states. The first column is the estimate from 
the Buchmiiller-Tye potential model as given in Ref. [11]. The second and third columns are the 
phenomenological results from electromagnetic annihilation decays for rric = 1.4 GeV at leading 
order (LO) and next-to-leading order (NLO) in ag- The errors are the statistical errors associated 
with the experimental inputs only. The bold-faced values are used in the predictions for the 
double-charmonium cross sections. 



potential model 

phenomenology 

H 


LO 

NLO 

Vc 

0.387 

0.222 ± 0.024 

0.297 ± 0.032 

Jl^ 

0.387 

0.208 ± 0.015 

0.335 ± 0.024 

r]c{2S) 

0.253 



i^{2S) 

0.253 

0.087 ± 0.006 

0.139 ± 0.010 

Xco(lP) 

0.107 

0.060 ± 0.015 

0.059 ± 0.015 

Xci(lP), hc{lP) 

0.107 



Xc2{lP) 

0.107 

0.033 ± 0.006 

0.053 ± 0.009 

V^i(iP) 

0.054 

0.095 ± 0.015 


V'2(1P), ??c2(lP) 

0.054 




III. NRQCD MATRIX ELEMENTS 

The ratios R for exclusive double-charmonium production calculated in Section II depend 
on the NRQCD matrix elements {Oi)h- In this section, we describe the phenomenological 
determination of these inputs. We also give estimates based on potential models. 


A. Potential models 


We can obtain estimates for the NRQCD matrix elements from the behavior of the 
wavefunctions near the origin in potential models. The expressions for the NRQCD matrix 
elements for P-wave, P-wave, and P-wave states are 


(Oi)s ~ 


(30a) 

(Oi)p 


(30b) 

{Oi)d - 


(30c) 


The values and derivatives of the radial wavefunctions at the origin for four potential models 
are given in Ref. [11]. Of these four potential models, the one that is most accurate at short 
distances is the Buchmiiller-Tye potential [12]. The values of the NRQCD matrix elements 
for this potential are given in the first column of Table I. 
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B. Phenomenology 


We can obtain phenomenological values for the NRQCD matrix elements and 

from the electronic decay rate of the J/ijj and from the photonic decay rate of the rj^ 
[1], The results for these decay rates, including the hrst QCD perturbative correction, are 


r[^c ^ 77] 
VlJ/i) e+e-] 


2e^7rQ;^ 

(Ql)j/V; A _ 8 ^ 

3 ml \ 3 TT / 


{Oi)vc A _ 20 -TT^ 
ml \ 6 T’’ / ’ 


(31a) 

(31b) 


We can obtain phenomenological values for the NRQCD matrix elements and (Oi )^^2 

from the photonic decay rates of the Xco and Xc 2 [!]• The results for these decay rates, 
including the hrst QCD perturbative correction, are 


r[xco - 

77] 

= Ge^TTo;' 

r[xc 2 - 

77] 

8 e^ 7 rQ; 

5 


(Qi)xco 

ml 

(Qi)xc 2 

ml 


( 371^-28 asV 

V 18 TT y ’ 

\ 2 
71 J 



(32a) 

(32b) 


The perturbative corrections in Eqs. (31a)-(32b) have been expressed as squares, because 
they can be calculated as corrections to the amplitudes. We can obtain a phenomenologi¬ 
cal value for the NRQCD matrix element {Oi)^ 2 iiD) from the electronic decay rate of the 
ipillD) = ^>(3770), which is in the same spin-symmetry multiplet as the ^>2(1-0); 


r[y>i(lD) —> e’''e 


18 m® 


(33) 


The QCD perturbative correction has not been calculated. 

We hrst determine the NRQCD matrix elements {Oi)h while neglecting QCD pertur¬ 
bative corrections. The experimental inputs are the electronic widths of the 7p{2S), 
and 7pi{lD), the photonic width of the rjc, and the widths and photonic branching fractions 
of the Xco and Xc 2 [13]. The only other input required is the charm quark mass me- The 
values of {Oi)h corresponding to me = 1.4 GeV are given in the column of Table I labelled 
LO. The error bars are the statistical errors associated with the experimental inputs only. 
To obtain {Oi)h for other values of me, we need to multiply the values in the Table by 
(mc/1.4 GeV)2+2^. 

We next determine the NRQCD matrix elements {Oi)h including the effects of QCD per¬ 
turbative corrections. We choose the QCD coupling constant to be as = 0.25 corresponding 
to a renormalization scale of 2me- The resulting values of {Oi)h for me = 1.4 GeV are given 
in the column of Table I labelled NLO. The error bars are the statistical errors associated 
with the experimental inputs only. To obtain {Oi)h for other values of me, we need to mul¬ 
tiply the values in the Table by (mc/1.4 GeV)^^^'^. Taking into account QCD perturbative 
corrections changes {Oi)h by a factor that ranges from 0.99 for Xco fo 1.61 for J/ip and 
Xc 2 - The NLO values are in closer agreement with the estimates from the Buchmiiller-Tye 
potential model than the LO values. 

One complication in the determination of the NRQCD matrix element for 7pi{lD) is that 
this state may have substantial mixings with the t/’(25') and also with continuum DD states. 
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If the mixing angle between 'ipi{lD) and is 0 and if mixing with continuum DD states 

is neglected, the expressions for the electronic decay rates of the and 'ipi{lD) are (31b) 

and (33) with the substitutions 


{Ol)^(2S) 


COS 0 (01)25^ ~ 

cos 0 (Oi)Jo + sin 


^(n 

Qnil 

VTE 


{Oi) 


ID 


1/2 

25 


(34a) 

(34b) 


A recent estimate of this effect suggests a mixing angle 0 = 12° [14], The resulting values 
of the NRQCD matrix elements are ( 01)25 = 0.095 GeV^ and {Oi)id = 0.013 GeV^. This 
value of (Oi)iD is about a factor of 7 smaller than the value of (Oi)^^(i£)) in Table I. Thus, 
if this mixing scenario is correct, the phenomenological estimate for {Oi)^ 2 {id) in Table I 
could overestimate cross sections for 02(10) by about a factor of 7. 

Within each spin-symmetry multiplet, the NRQGD matrix elements {Oi)h should have 
differences of order which we expect to be about 30%. We choose to use the most 
precise phenomenological value within each spin-symmetry multiplet for all members of 
that multiplet. Specihcally, we use the bold-faced values in the NLO column of Table I for 
S'-wave and P-wave states and we use the bold-faced value in the LO column for the P-wave 
states ? 7 c 2 ( 1 P), 0 i(1 P), and 02 ( 1 P). 


IV. RELATIVISTIC CORRECTIONS 

In this section, we calculate the relativistic corrections to the cross sections for the S'-wave 
double charmonium. We also give a phenomenological determination of the NRQCD factors 
{v‘^)h that appear in those relativistic corrections. 


A. NRQCD factor 

The leading relativistic correction to the GSM amplitudes for a charmonium H are con¬ 
veniently expressed in terms of a quantity that is denoted by {v‘^)h in Ref. [1]. It can be 
dehned formally as a ratio of matrix elements in NRQCD. For example, in the case of r]c, it 
can be written as 


(0|xt(-iD)2^|^0 


(35) 


The naive interpretation of {v‘^)h is the average value of for the charm quark in the 
charmonium H. In the Buchmiiller-Tye potential model, the average value of weighted by 
the probability density is 0.23 for the IS* states rjc and J/0 and 0.29 for the 2S states rfci^S) 
and 0(25*) [12]. However the proper interpretation of the ratio of matrix elements in (35) is 
the average value of weighted by the wavefunction. Unfortunately this quantity has power 
ultraviolet divergences and requires a subtraction. For example, the wavefunction for the 
IS states in potential models can be approximated fairly accurately by a momentum space 
wavefunction of the form 0 (p) = l/(p^-|-m^uf^)^ where U 15 is a phenomenological parameter. 
The integral J d^pp^'ip{p) has a linear ultraviolet divergence. Minimal subtraction of this 
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linear divergence gives the negative value {v'^)h = Thus the extraction of estimates 

of {v‘^)h from potential models is not straightforward. 

There is a connection between the quantity {v‘^)h and the mass of the charmonium state 
H that was first derived by Gremm and Kapustin [15]. The most convenient form of this 
relation for relativistic applications is 


Ml = Ami (l + + ■■■), (36) 

where the corrections are of order mlv*. The mass me that appears in this relation is the 
pole mass. The pole mass suffers from renormalon ambiguities, but those ambiguities are 
largely compensated by corresponding ambiguities in the matrix elements that dehne {v‘^)h 
[16, 17]. 

We can use the Gremm-Kapustin relation (36) to obtain a phenomenological determina¬ 
tion of {v'^)h using the mass Mh of the charmonium state as input; 



Ml - Ami 

Ami 


(37) 


where me is the pole mass of the charm quark. The masses for the charmonium states rje, 
J/iIj, and are well-measured. The r]e{2S) was only recently discovered by the BELLE 

Gollaboration with a mass of M^^(^ 2 S) = 3654 ± 6 ± 8 MeV [18]. If we set me = 1.4 GeV, the 
values of {v‘^)h for the S'-wave states are {v‘^)nc = 0.13, {v‘^)j/^ = 0.22, {v'^)rj^( 2 S) = 0.70, and 
{v‘^)tp{ 2 S) = 0.73. The values of {v‘^)h for the 2S states are uncomfortably large, but those 
large values are necessary to compensate for the fact that 2me = 2.8 GeV is far from the 
mass of the 2S states. 


B. Relativistic correction factor 


The relativistic correction to the cross section for the process J/ip + rje can be calculated 
by replacing the amplitude factors (9a) and (13a) by 


M[t]c] 




Vc 


ANeml{l + (u^) 


VcJ 


1/2 


m: 


A + Car 1° 




(38) 

(39) 


The prefactors take into account the relativistic normalizations of cc states and charmonium 
states. Strictly speaking, the factor of (1 -|- should be expanded out to first 

order in {v‘^)h- However if we use phenomenological determinations of the NRQGD matrix 
elements {Oi)h, the prefactor in (38) or (39) cancels. We therefore choose not to expand 
the prefactors. 

There are relativistic corrections to the electromagnetic annihilation decay rates used to 
determine the NRQGD matrix elements in Table 1. For the decay rates of the r]e and the 
J/ip given in Eqs. (31a) and (31b), the leading relativistic correction can be expressed as a 
multiplicative factor 


1 


1 2\ \ ^ Mh 2me 

6 ^ ^ 2me{l + {v'^) h) ^ Mh' 


(40) 


13 



The correction has been expressed as the product of three factors. The first factor, which 
appears squared, conies from the expansion of the amplitude in powers of the relative velocity 
of the cc pair. The second factor comes from the prefactor in Eq. (38) or (39). The last 
factor comes from the relativistic normalization factor 1/{2Mh) in the standard expression 
for the decay rate. Note that the factors of Mh cancel in Eq. (40). 

We have calculated the relativistic corrections to the cross section for J/^j + fjc- The 
leading correction to the ratio R in Eq. (18) can be expressed as the multiplicative factor 


1 + 


8F + 3(F + 4)r2 - 2 \ 2Y + (Y + Uy - 5r 2 


12(r2 - Y) 


-{v )j/^ + 


i2y-Y) 


v^) 


Vc 


M 


X 


j/i/. 


M, 


nc 


2mc{l + y)j/^) 2mc{l + y)riJ 


X 


1 3/2 


2(1 -r2) 




. (41) 


The correction has been written as the product of three factors. The first factor, which 
appears squared, comes from the expansion of the amplitude in powers of the relative velocity 
of the cc pair. The second factor comes from the prefactors in Eqs. (38) and (39). The last 
factor is the nonrelativistic expansion of the term (PoMZ-Sbeam)^ divided by its value in the 
nonrelativistic limit, where one power is from the phase space factor (20) and the other two 
are from the square of the amplitude in Eq. (16). 

Our final result for the relativistic correction can be expressed as a multiplicative factor 
obtained by dividing (41) by a factor of (40) for each of the charmonium states. We express 
it in the form: 



8 F + 3(F + 4y - 5r^ 
12(r2 -F) 




X 




2Y + {Y + uy - 5r^ 

12(r2 - Y) 

f Pcm/ -^beam \ ^ 
4 m 2 V (1 ~ 7-2^ 1/2 j 



(42) 


V. PREDICTIONS FOR B FACTORIES 

In this section, we calculate the cross sections for exclusive double-charmonium produc¬ 
tion in e’''e“ annihilation at the B factories. We also give a careful analysis of the errors in 
the predictions for J/'ip + rjc- 


A. Cross sections 

The results in Section II were expressed in terms of the ratio R dehned in Eq. (2). The 
corresponding cross sections are 

a[H, + H 2 ] = ^R[Hi + H 2 ]. (43) 

3s 

The ratios R depend on a number of inputs: the coupling constants and a, the charm 
quark mass rric, and the NRQCD matrix elements {Oi)h- 
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TABLE II: Cross sections in fb for e+e annihilation into double-charmonium states Hi + H 2 
without relativistic corrections. The errors are only those from variations in the NLO pole mass 
rric = 1.4 ± 0.2 GeV. 


H 2 \ Hi 

Jjlp 

iP{2S) 

hc{lP) 

iPi{lD) 

iP2{ID) 

Vc 

3.78 

± 

1.26 

1.57 

± 

0.52 

0.308 

± 

0.017 

0.106 ± 0.025 

1.04 ± 

0.23 

r}c{2S) 

1.57 

± 

0.52 

0.65 

± 

0.22 

0.128 

± 

0.007 

0.044 ± O.OlO 

0.43 ± 

0.09 

Xco(l-P) 

2.40 

± 

1.02 

1.00 

± 

0.42 

0.053 

± 

0.019 




Xci(lF) 

0.38 

± 

0.12 

0.16 

± 

0.05 

0.258 

± 

0.064 




Xc2{lP) 

0.69 

± 

0.13 

0.29 

± 

0.06 

0.017 

± 

0.002 




??c2(1D) 

0.35 

± 

0.05 

0.14 

± 

0.02 








The value of the QCD coupling constant as depends on the choice of the scale /i. In the 
QCD diagrams of Fig. 1, the invariant mass of the gluon is ^/i/2. We therefore choose the 
scale to be /i = 5.3 GeV. The resulting value of the QCD coupling constant is as{fi) = 0.21. 

The numerical value for the pole mass rric of the charm quark is unstable under perturba¬ 
tive corrections, so it must be treated with care. Since the expressions for the electromagnetic 
annihilation decay rates in Eqs. (31a)-(32b) include the perturbative correction of order o;*, 
the appropriate choice for the charm quark mass rric in these expressions is the pole mass 
with corrections of order ag included. It can be expressed as 


rric 


mc{mc) 



(44) 


Taking the running mass of the charm quark to be rfic{ffic) = 1.2 ± 0.2 GeV, the NLO pole 
mass is rric = 1-4 ± 0.2 GeV. 

Our predictions for the double-charmonium cross sections without relativistic corrections 
are given in Table II. The error bars are those associated with the uncertainty in the NLO 
pole mass rric only. 

Our predictions for the double-charmonium cross sections for the S'-wave states including 
the leading relativistic correction are obtained by multiplying the values in Table II by the 
factor (42). We use the values of {v‘^)h obtained from Eq. (37), which follows from the 
Gremm-Kapustin relation. The resulting cross sections are given in Table III. The error 
bars are those associated with the uncertainty in the NLO pole mass rric only. The relativistic 
corrections increase the central values of the cross sections by about 2 for J/ilr + rjc, by about 
5 for J /i/j + ric{2S) , by about 4 for +ric, and by about 8 for iIj{2S) +ric{2S). Although 

the total correction factor for J/'ip+r]c is signihcantly larger than 1, it is the product of several 
modest correction factors that all go in the same direction. The largest individual correction 
factor for J/ip+rjc is (1.28)^ coming from the expansion of the amplitnde. The corresponding 
factors for J / ip + r]c{2S)^ i’{2S) + r]c^ and ip{2S) + r]c{2S) are (1.80)^, (1.64)^, and (2.16)^, 
respectively. These large correction factors indicate that the relativistic corrections to the 
cross sections involving 2S states are too large to be calculated reliably using the method 
we have chosen. 

Note that onr method for calculating the relativistic correction signihcantly increases the 
sensitivity to the charm qnark mass. The errors from varying rric in Table II are about 50% 
for the S'-wave states, while the errors in Table III correspond to increasing or decreasing 
the cross section by abont a factor of 3. The strong sensitivity to rric is another indication 
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TABLE III: Cross sections in fb for e'’'e“ annihilation into 5-wave double-charmonium states 
Hi + H 2 including relativistic corrections. The errors are only those from variations in the NLO 
pole mass rric = 1.4 ± 0.2 GeV. 


H2 \ Hi 

Jli, 

i/.(25) 

Vc 

f-7 /j+10.9 
' *^-4.1 

Fill! 

r]c{2S) 

'•'3-4.1 



that our method for calculating the relativistic corrections is unreliable. We will therefore 
take the values in Table II to be our predictions for the cross sections and use Table III as 
an indication of the possible size of the relativistic corrections. 


B. Perturbative corrections 

The QCD perturbative corrections to the electromagnetic annihilation decay rates used 
to determine the NRQCD matrix elements in Table I have already been taken into account. 
The QCD perturbative corrections to the cross section for J/ip + rjc have not yet been calcu¬ 
lated. However parts of the perturbative corrections are related to perturbative corrections 
that have been calculated and we can use these to estimate the order of magnitude of the 
perturbative corrections. 

Some of the perturbative corrections can be associated with the wavefunctions of the 
Tjc and J/'ip. For the QED diagrams in Fig. 2, the QCD perturbative corrections would be 
very closely related to those in the electromagnetic annihilation decay rates (31a) and (31b). 
However for the QCD diagrams in Fig. 1, the QCD perturbative corrections associated with 
the wavefunction could be very different. In the expressions for electromagnetic annihilation 
decay rates in Eqs. (31a)-(32b), we have 4 examples of perturbative corrections associated 
with wavefunctions. The root-mean-square of the 4 coefficients of is 0.66. We will 
therefore take (1 ± 0.660:5)^ as our estimate for the perturbative correction associated with 
each charmoninm wavefunction. 

There are other pertnrbative corrections that can be associated with the electromagnetic 
charm current cy^c. As an estimate for the magnitnde of these corrections, we can nse the 
perturbative correction to the inclusive charm cross section. The corresponding ratio R is 
[19] 


R[cc -I- X] = 3e^(l -h r^/8)(l — 


, , . dr"' \ as 

4 16 / TT 


(45) 


If we set nic = 1.4 GeV and Ebeam = 5.3 GEV, then = 0.28 and the perturbative correction 
gives a mnltiplicative factor (1 -|- O.lOo;.,)^. 

There are also perturbative corrections that can be associated with the QGD coupling 
constants. We can estimate the size of these pertnrbative corrections by varying the scale /i 
up or down by a factor of 2. The factor in the cross section changes by a mnltiplicative 
factor of (1 ± 0.920:5)^. 

To obtain an estimate of the errors in the double charmoninm cross section associated with 
perturbative corrections, we will add in qnadrature the coefficients of a* in the perturbative 
correction factors associated with each wavefnnction, the charm current, and the factors of 
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as- The resulting correction is a multiplicative factor (1 ± If we set as = 0.21, this 

factor ranges from about 0.53 to about 1.62. Thus we should not be surprised if the QCD 
perturbative corrections changed the predictions by 60%. 


C. Color-octet contributions 


According to the NRQCD factorization formalism [1], the inclusive double-charmonium 
cross sections are obtained by replacing the decay NRQCD matrix elements in the color- 
singlet model terms by production NRQCD matrix elements and by adding additional terms 
involving color-octet matrix elements. The color-octet contributions at this order in ag can 
be obtained from the results in Section II by replacing the NRQCD matrix elements by 
appropriate color-octet matrix elements and by replacing the constants X and Y by the 
constants Xg = and Ig = —3. 

As an illustration, we consider the inclusive production of J/'ip + fjc- The leading color- 
octet contribution for J/ip + 7]^ + X can be obtained from the color-singlet model result for 
Xcj + he in Eq. (25) by substituting {Oi{^Pj))^^j {0^^'^{^Pj))/{2J + 1), {Oi{^Pi))h, 

(Og“(^Pi)), X ^ Xg and Y Yg. The ratio R for J/ip + rie + X at this order includes two 
terms: 


R[J/^ + Vc + X] 




9 


3ml 


m-i 


(46) 


By the velocity-scaling rules of NRQCD [1], the color-octet term is suppressed by a power 
of n®, but that is partly compensated by an enhancement factor of 1/r®. Omitting the 
color-octet term and applying the vacuum-saturation approximation to the NRQCD matrix 
elements in the color-singlet term, we recover the nonrelativistic limit of the ratio R in 
Eq. (18) for the exclusive J/ip + rje hnal state. 

Color-octet processes can also contribute to the exclusive cross section for J/'ip + rjc. The 
2 gluons emitted by one color-octet cc pair in the transition to a color-singlet state that 
can form charmonium can be absorbed by the other color-octet cc pair. The amplitude is 
suppressed by relative to the color-singlet amplitude. The leading contribution to the 
cross section will come from the interference between these two amplitudes and so will be 
suppressed only by The interference terms cancel when summed over all possible hnal 
states, so they do not appear in the inclusive cross section (46). The color-octet contributions 
will also have suppression factors of that guarantee consistency with the helicity selection 
rules of perturbative QCD in the limit r —0. 


D. Phenomenology 

The BELLE Collaboration has recently measured the cross section for J/ip + rjc [4]. The 
J/ip was detected through its decays into and e+e“, which have a combined branching 

fraction of about 12%. The pc was observed as a peak in the momentum spectrum of the 
J/ip corresponding to the 2-body process J/ip + rjc. The measured cross section is 

a[J/iP + 7 ]e] X B[> 4] = {33PI ± 9) fb, (47) 
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where B[> 4] is the branching fraction for the rjc to decay into at least 4 charged particles. 
Since B[> 4] < 1, the right side of Eq. (47) is a lower bound on the cross section for J/tjj + ric. 

The lower bound provided by Eq. (47) is about an order of magnitude larger than the 
central value 3.8 fb of the calculated cross section for J/ip + rjc in Table II. The largest 
theoretical errors are QCD perturbative corrections, which we estimate to give an uncertainty 
of roughly 60%, the value of nic, which we estimate to give an uncertainty of roughly 50%, 
and a relativistic correction that we have not been able to quantify with conhdence. If 
we take the calculations of the relativistic corrections in Table III seriously, the extreme 
upper end of the prediction is marginally compatible with the BELLE measurement. In our 
further discussion, we will ignore the large discrepancy between the predicted cross section 
for J/tjj + rjc and the BELLE measurement. We will focus on predictions for the ratios of 
cross sections from Table II under the assumption that many of the theoretical errors will 
cancel in the ratio. 

In addition to measuring the cross section for J/'^ + rjc, the BELLE Collaboration saw 
evidence for J/ip + rjci^S) and J/ip + Xcoi^P) [4] events. A 3-peak £t to the momentum 
spectrum of the J/'^ gave approximately 67, 42, and 39 events for rjc, ric{2S), and Xcoi^P) 
with fluctuations of 12-15 events. The proportion of J/ip -f rjci J/p) -7 ric{2S), and J/ip + 
Xcoip-P) events, 1.00 : 0.63 ± 0.25 : 0.58 ± 0.24, is consistent with the proportions 1.00 : 
0.41 : 0.63 of the cross sections in Table 11. These proportions are insensitive to the choice 
of rric. The absence of peaks corresponding to Xai/^P) and Xc 2 {XP) is also consistent with 
Table 11. The cross sections for J/ip + Xcj{^P) for J = ^ and 2 are predicted to be smaller 
than for J = 0 by factors of about 0.16 and 0.29, respectively. 

If the cross sections for the narrow H-wave states are large enough, they could be discov¬ 
ered at the B factories. The state r]c2{^D) could be observed as a peak in the momentum 
spectrum oi J/ip corresponding to the 2-body process J/ip + ?7c2(171). The prediction in 
Table II for the cross section for J/ip + ?7c2(171) is smaller than that for J/ip + rjc by about 
a factor 0.09. It might also be possible to discover the H-wave state '^2(1^*) as a peak in 
the momentum spectrum of r]c corresponding to the 2-body process ip 2 {lD) + rjc- The rjc 
could be detected through its decay into KKn, whose branching fraction is about 6%. The 
prediction in Table II for the cross section for ip 2 {lD) +r]c is smaller than that for J/ip + t]c 
only by about a factor 0.27. Our predictions for the cross sections for J/ip -7 ?7c2(171) and 
ip 2 p\-D) -7 rjc are based on a phenomenological determination of the NRQCD matrix ele¬ 
ments that ignored mixing between the ip{2S) and the ipi{lD). If there is significant mixing 
between these two states, the cross sections for S'-wave -7 H-wave could be a factor of 7 
smaller. 

In summary, we have calculated the cross sections for e’''e“ annihilation into exclusive 
double charmonium states with opposite charge conjugation. Many of the cross sections are 
large enough to be observed at B factories. In particular, it may be possible to discover the 
D-wave states r]c 2 {lD) and ip 2 {lD). The cross sections for double charmonium suffer from 
fewer theoretical uncertainties than inclusive charmonium cross sections. The largest uncer¬ 
tainty comes from relativistic corrections. Measurements of exclusive double charmonium 
cross sections will provide strong motivation for developing reliable methods for calculating 
the relativistic corrections to quarkonium cross sections. 

Note added: Liu, He, and Chao have calculated the a‘^a‘^ terms in the cross sections 
for annihilation into J/ip + H, H = ric,Xco,Xci,Xc2 [20]. Their results are consistent 
with ours. In Ref. [21], Liu, He, and Chao calculated the a'^a^ terms in the differential cross 
sections for J/ip + rjc, J/ip + Xco, and rjc -7 he- Their results for J/ip + r]c and J/ip + Xco agree 
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with ours. Our corrected result for r]c + he in Eq. (23) agrees with their result. 


APPENDIX A; ANGULAR DISTRIBUTIONS 

In this appendix, we give the angular distributions dR/dx for iJi(Ai) + H 2 {\- 2 ) 

for each of the helicity states that contribute at order a‘^a‘1 or The angular variable is 
X = cos 6, where 6 is the angle between e~ and Hi in the center-of-mass frame. The 
results for R in the text are obtained by summing over all the helicities and integrating over 
-1 < a; < +1. 


1. S-wave + S'-wave 

The angular distribution for J/ip + rjc is 

dR ^ 0\Q/o {.hdl) JJ^i^Ol) r]e 


^|J/V(±l) + dJ = ^ XHr^- - Y)V(1 - 


rrvi 


(l + x^), (Al) 


The cross section for the longitudinal helicity component Ai = 0 of the J/ijj vanishes. 


2. S-wave -|- P-'wseve 

The angular distributions for J/'ip + Xcj are 

:V(i - a2,x). 


flW(A.)Vw(V)| = ^V 




The non-vanishing entries of Fj{Xi, X 2 ,x) are 


(A2) 


Fo(0,0,a;) = ^r^[4 + 2(F + 5)r^ — 3r"^]^(l — 

(A3a) 

Fo(±l,0,a;) = l[AY - Q{Y + + x^), 

O 

(A3b) 

Fi(0,±l,a;) = ^r^[2(3F + 4)-7r2]2(l + a;2), 
lb 

(A3c) 

Fi(±l,0,a;) = ^[SY-2Yr^+ r^f{l + x^), 
lb 

(A3d) 

Fi(±l,±l,a;) = \r‘^[AY+ 2{Y+ 2)r‘^-?>r^f{l-x^), 

8 

(A3e) 

F2(0,0,a;) = ^r^[4 + 2(F — l)r^ — 3r^]^(l — 

(A3f) 

F2(0,±l,a;) = ^r"[2(F + 2)-5r2]2(l + a;2), 
lb 

(A3g) 

F2(±l,0,a;) = 4[8y-b(y+ 2)r2 + llr^]2(l + a;2), 
lb 

(A3h) 

F2(±l,±l,a;) = ^r2[4F-2(F + 2)r2 + 3r^]^(l-a;^), 

8 

(A3i) 

F 2 (±l,± 2 ,a;) = lr^[2Y - {1 + x^). 

O 

(A3j) 
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The angular distribution for rjc + he is 


where 


+ 144 




V(1 - r")‘''"ihh!4hhi//(A, 

m° 


X) 


hf(0,a;) = ^(1 — a;^)(3r‘^ — 6r^ + 4)^, 

H{±.l^x) = -(1 + a;^)r^(r^ — 2)^. 

8 


(A4) 


(A5a) 

(A5b) 


3. P-wave + P-wave 


The angular distributions for he + Xcj are 

+ X.i(A 2)] = Alk<2lkiG,(A,, A2 ,i). 

ax 108 mix 


(A6) 


where non-vanishing entries of Gj(Ai, \ 2 -,x) are given by 

Go(±l, 0 ,a;) = {% - r'^f {I + x^), 

O 

G'i(0,0,a;) = 18(1 

Gi(0,±l,a;) = ^r^(2-r^)^(l+ a;^), 
lb 

Gi(±l, 0 ,a;) = -^r^( 2 -r^)^(l-f-a;^), 
lb 

6 * 2 ( 0 , ±l,a;) = -^r^(4-5r^)^(l-f 
lb 

G'2(±l,0,a;) = ^r®(l + a;^), 
lb 

G2(±l,±2,a;) = ^r®(l-Fa;^). 


(A7a) 

(A7b) 

(A7c) 

(A7d) 

(A7e) 

(A7f) 

(A7g) 


Note that G2(±l,0,a;) is more suppressed than the prediction from the hadron helicity 
conservation rule. 


4. S-wave -|- P-wave 

The angular distributions for J/ip + rie 2 are 

^|J/V.(±l) + 1,02(0)1 = ^ X^Y - 2rrrHl - (1 + (A8) 

The cross sections for the longitudinal helicity component Ai = 0 of the J/ijj and the helicity 
components A2 = ±1 and ±2 of the ?7c2 vanish. 
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The angular distributions for + rjc are 






X 2 ( 26 y - 21r2 + lOr^r^l - + a;2 


^).(A9) 


The cross section for the longitudinal helicity component Ai = 0 of the 'ipi vanishes. 
The angular distributions for ip 2 + Vc are 


f [V-.(A.) + »ol = ^ A'V(1 - 


X) 


(AlO) 


The non-vanishing entries of H{Xi,x) are 


9 


H{0,x) = ), 

Hi±l^x) = -^r^(7 — 4r^)^(l-I- 
16 


(Alla) 

(Allb) 
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Erratum: Exclusive double-charmonium production 
from e+e“ annihilation into a virtual photon 
[Phys. Rev. D 67, 054007 (2003)] 
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In Ref. P|, we calculated the cross sections for exclusive double-charmonium production from e^e“ annihilation 
into a virtual photon. We made an error in the calculation by omitting a relative minus sign between diagrams that 
differ by the interchange of identical fermions in the final state. That error was also made in Refs. 01 and 0. In 
Ref. 0, the error was the omission of a relative minus sign between the sum of the four Feynman diagrams in Fig. 1 
and the sum of the two QED Feynman diagrams in Fig. 2. We also made a separate error in the calculation of the 
cross section for ijc he- We show corrections in boldface type, except in display equations. 

The expression in Eq. (17) for the coefficient A in the amplitude for J/i/; -|- rjc should be corrected as follows: 


A = 


1287ra= 


{{Ol) J/'4i{Ol)r)„) 


1/2 


N^-l 

2m 


, 1 ec«'\ 

NcUs r2 as ) 


(17) 


The sign in front of 1/r^ term has been reversed. In the expression for the coefficient Y in Eq. (19b), the overall sign 
should also be reversed: 


Y = - 


a 


OLs 



(19b) 


In the subsequent sentence, the numerical value of Y should have the opposite sign: Y = —0.0344. 

In Section II, the descriptions for the changes in various cross sections due to the QED contribution should be 
corrected as follows: 

• Sec. II C (S-wave -I- S-wave) above Eq. (20): If we set ^/s = 10.6 GeV and rric = 1.4 GeV, the electromagnetic 
correction increases the cross section by 29%. 

• Sec. II D (S-wave -I- P-wave) above Eq. (23): If we set -y/s = 10.6 GeV and TOc = 1-4 GeV, the QED corrections 
change the cross sections by -|-5.0%, —5.5%, and +11% for J = 0, 1, and 2, respectively. 

• Sec. II F (S-wave + D-wave) below Eq. (26): If we set ^/s = 10.6 GeV and rric = 1.4 GeV, the QED correction 
increases the cross section by about 15%. 

• Sec. II F (S-wave + D-wave) below Eq. (27): The QED contribution increases the cross section by 41%. 

There are several corrections to the text in Section V: 


• The relativistic corrections increase the central values of the cross sections by about 2 for J/ip -\- ric, by about 
5 for J/tp + ric{2S), by about 4 for 'ip{2S) + rje, and by about 8 for ip{2S) + r]c{2S). 

• The largest individual correction factor for J/ip ijc is (1.28)^ coming from the expansion of the amplitude. 
The corresponding factors for J/pj -\- r]c{2S), tp{2S) + +, and ipiflS) -\-r]c{2S) are (1.80)^, (1.64)^, and (2.16)^, 
respectively. 

• The lower bound provided by Eq. (46) is about an order of magnitude larger than the central value 3.8 fb of 
the calculated cross section for J/ip r]c in Table II. 

• The proportion of [J/ip + r]c, J/ip + +(2S'), and J/ip + Xco(l^’) events,] 1.00 : 0.63 ± 0.25 : 0.58 ± 0.24, is 
consistent with the proportions 1.00 : 0.41 : 0.63 of the cross sections in Table im 

• The cross sections for J/ip + Xcj{^P) for J = 1 and 2 are predicted to be smaller than for J = 0 by factors of 
about 0.16 and 0.29, respectively. 
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TABLE II: Cross sections in fb for e^e annihilation into double-charmonium states Hi + H 2 without relativistic corrections. 
The errors are only those from variations in the NLO pole mass rric = 1.4 ± 0.2 GeV. 


H 2 \ Hi 

J/^ 

V>(25) 

heilP) 

V'i(lD) 

M^D) 

Vc 

3.78 ± 1.26 

1.57 ± 0.52 

0.308 ± 0.017 

0.106 ± 0.025 

1.04 ± 0.23 

vy2S) 

1.57 ± 0.52 

0.65 ± 0.22 

0.128 ± 0.007 

0.044 ± 0.010 

0.43 ± 0.09 

Xco(lP) 

2.40 ± 1.02 

1.00 ± 0.42 

0.053 ± 0.019 



Xci(lR) 

0.38 ± 0.12 

0.16 ± 0.05 

0.258 ± 0.064 



Xc2(lP) 

0.69 ± 0.13 

0.29 ± 0.06 

0.017 ± 0.002 



Vc2{1D) 

0.35 ± 0.05 

0.14 ± 0.02 





TABLE III: Cross sections in fb for e"'"e annihilation into S-wave double-charmonium states Hi -|- H 2 including relativistic 
corrections. The errors are only those from variations in the NLO pole mass rric = 1.4 ± 0.2 GeV. 


H 2 \ Hi 

J/V> 

V>(25) 

Vc 

7 4+10-9 
* *^-4.1 

Q 4 + 9-5 

Vc{2S) 

7 g + 11-8 
* *1^-4.1 

K 0 + 9.1 
^•^-2.9 


• The prediction in Table HII for the cross section for J/V' + r]c 2 {lD) is smaller than that for J/'0 + t]c by about a 
factor 0.09. 


• The prediction in Table ITU for the cross section for 'ip 2 {^D) + rjc is smaller than that for J/ip + rjc only by about 
a factor 0.27. 

In Ref. P|, we also made an error in the calculation of the cross section for r]c + he- The result that was obtained 
for the ratio R in Eq. (23) was zero. The correct result is 


where the function H(r) is 




144 


H{r) = 2r^{r^-2y + {3r*-6r^+4y. 


(23) 


The dependence on a appears only in the overall factor because the QED contribution comes from diagrams with 
the same topology as the (^D diagrams in Fig. 1 of Ref. 0. The QED contribution from the photon-fragmentation 
diagrams in Fig. 2 of Ref. [l| vanishes because the -I- parity of he does not allow the direct coupling to a single virtual 
photon. The pure QCD contribution to the result in Eq. lP|) is obtained by substituting X = 4/9. This result was 
first calculated correctly in Ref. Q . The angular distribution for rjc + he is 


where 


dR 

dx 


\i. + (..(A)! = 


{Ol)riAOl)h^ 


144 


7L(A, x) 


H{0,x) = ^(l-x2)(3r‘‘-6r2-b4)2, 

H{±i,x) = lii + xyr^y^ -2y. 
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The error in the calculations produced errors in the numerical predictions for some of the cross sections given in 
Tables II and III of Ref. P|. Those Table are reproduced here, with the corrected values given in boldface. 

We thank Chaehyun Yu for his help in checking some of the calculations. 
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